CERN-TH/97-240 



|hep-ph/9709327 



Theoretical Analysis oi B — > D**7t Decays 



M. Neubert 

Theory Division, CERN, CH-1211 Geneva 23, Switzerland 



Abstract 

The decays B D**tt, where D** is a narrow p-wave charm resonance, are in- 
vestigated in the context of a generahzed factorization hypothesis, including the 
leading nonfactorizable corrections. The decay rates for these processes are related 
to the corresponding semileptonic rates at maximum recoil. It is pointed out that 
the current data on the branching ratio for the decay vr^ may pose a 

problem for theory. We predict that future, more accurate measurements of this 
branching ratio will find a value B(i?~ — > D^^^ -k~) 4 x 10~^, a factor 5 lower 
than the current central value. 
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1 Introduction 



Recently, the CLEO Collaboration has reported the first measurement of the branching 
ratios for the decays B~ D^7r~ and B~ -k^ , where Di and D2 (generically 

denoted as D**) are p-wave charm mesons. In the heavy-quark limit, these are the 
members of a spin-symmetry doublet in which the light degrees of freedom have total 
angular momentum J = §• The measured branching ratios are 

B(5^ ^ TT") X B(L)° ^ 0*+-^-) = (7.8 ± 1.9) X 10^^ , 
B{B- Df -n-) X B{Df D*+n-) = (4.2 ± 1.7) x 10^^ . (1) 

The purpose of this letter is to interpret these observations together with the existing data 
on the semileptonic decays B~ —>■ D^i~u and B~ D2^i~u. For the corresponding 
branching ratios, the CLEO Collaboration quotes [0| 

B{B- D° ru) X B(D° ^ D*+7i-) = (0.37 ± 0.10)% , 
B{B- Dfrv) X B{Df D*+n-) = (0.059 ± 0.067)% 

< 0.16% (90% CL) . (2) 
These numbers are consistent with results reported by the ALEPH Collaboration 
B{B Di i-jjX) X B(Di ^ D*TT^) = (0.49 ± 0.11)% , 
B{B DfruX) X B(D2*° ^ D*+7i-) < 0.34% (95% CL) , 
B{B DfruX) X B{Df < 0.33% (95% CL) , 

B{B D*+ ruX) X B{Dl+ 7r+) < 0.26% (95% CL) , (3) 

if one assumes that in most cases no additional particles X are produced. Semileptonic 
B decays into p-wave charm mesons are important with regard to understanding how 
the total inclusive semileptonic branching ratio of B mesons is composed out of exclusive 
modes. The ALEPH Collaboration finds a sizable branching ratio of (2.26 ± 0.44)% for 
the sum of all (resonant or nonresonant) decays B D^*^7ii~h' [Q]. These channels 
make up for about 20% of the total semileptonic branching ratio. 

Unfortunately, the D** D^*^7r branching ratios entering the above measurements 
have not yet been determined experimentally. Under the assumption that all Di mesons 
decay into D*n (decays into Dn are forbidden by angular momentum conservation), 
and using isospin invariance, one obtains B{D^ —>■ D*"'"7r~) = 2/3. For D2 mesons, 
the ratio B(D;° D+n~)/B{Df D*+tt-) = 2.3 ± 0.8 has been measured @, |]. 
Under the assumption that D2 — > D tt and D2 -D*7r are the only decay modes, 
and using isospin invariance, one then obtains B{D2^ — > D*"'"7r~) = (20 ± 5)% and 
B{Df D+TT-) = B{D*2^ 7r+) = (47 ± 5)%. 

In the case of neutral B mesons, there is a simple relation between hadronic decays 
involving a negatively charged pion in the final state and the corresponding semileptonic 
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decays with the pion replaced by a lepton-neutrino pair 0. In general, the amplitude 
for the hadronic decay H^tt^ , where represents an arbitrary charm meson (or 

a sum of hadrons with total charm number one), can be written as 

^had = ^ V,bV:,a, (vr-l J7m(1 - Ts)^ I ) {H^\ C7^(l - 75)^ \B') , (4) 

where oi parametrizes the complicated hadronic matrix elements of four-quark operators, 
which enter the theoretical description of nonleptonic amplitudes. This will be discussed 
in more detail below. The squared amplitude takes the form 

lAadP = ^ \VctV:f \a,\'flq,q.H^\ (5) 
where q = p — p' is the pion momentum, f^^ the pion decay constant, and 

H^^- = {H+\ era - 75)^ \B') {B^\ lYil - 75)c m) (6) 

the "hadronic tensor" . The same tensor appears in the description of semileptonic decays. 
Indeed, the squared amplitude for the process i'u is given by 

\A,,\' = ^\V^\'L,,H^\ (7) 

where L^,^ is the leptonic tensor. In the special case where = {pe + p^Y = 0, the 
leptonic tensor takes the simple form 



16x(l - x)qf,qu , (8) 



where x = Ei/ Ef^^^ is the scaled lepton energy in the rest frame of the B meson. Thus, 
in the limit where the pion mass is neglected, the nonleptonic decay rate is proportional 
to the semileptonic rate at q^ = 0. Working out the trivial phase-space factors, we obtain 
Bjorken's relation 

- = Qnf^\Vud\ \ai\ +0\—j\ (9) 



dr(50 ^ H+i'u)/dq^ 



Ml, 



irrespective of the nature of the hadron state He. The corrections of order M^/M^ have 
been worked out for the cases He = D and D* and are found to be negligible §]. Relation 
(y) still contains an unknown hadronic parameter Oi. However, for energetic two-body 
decays such as B^ D^*^^ii^ one can argue that the value of this parameter is close to 
unity. The reason is that the pion has a large energy in the rest frame of the decaying 
B meson. Once its constituent quarks have grouped together in a colour-singlet state, 
they form a fast-moving colour dipole which decouples from long-wavelength gluons. 
Only hard gluons with virtualities of order Mb are effective in rearranging the quarks. 
This "colour transparency argument" |^ suggests that the nonfactorizable contributions 
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to the hadronic decay amplitude are switched off at low momentum scales. Hence, one 
expects that ai ~ 1 for energetic two-body decays. This assertion is further strengthened 
by the 1/Nc expansion, which shows that ai = 1 + 0{1/N^) [Q. A determination of ai 
from B'^ D^*^^TT~ decays gives |ai| = 1.08±0.11, supporting the theoretical arguments 
just presented. A similar value is expected to apply for two-body decays into p-wave 
charm mesons. 

Unfortunately, the situation is more complicated for the decays of charged B mesons. 
Whereas the semileptonic decay rates of B~ and B^ mesons are related to each other by 
isospin invariance, this is not so for the hadronic rates. We will now discuss how relation 
(D must be modified in this case. 

2 Hadronic decay amplitudes for B — > D**7T 

The part of the effective weak Hamiltonian relevant to b cud transitions is given by 
Hefi = ^ V^tV:, {ci(/x) idu)icb) + cM {cu){dh) + ...}, (10) 

where {du) = d'y^{l — 75)^ etc. are left-handed, colour-singlet quark currents. The 
Wilson coefficients Cj(/i) are known to next-to-leading order. At the scale fi = rrih, they 
have the values ci(mfe) ^1.1 and C2(mb) ^ —0.3. These coefficients take into account 
the short- distance corrections arising from the exchange of hard gluons. The effects of 
soft gluons (with virtualities below the scale /i) remain in the hadronic matrix elements 
of the local four-quark operators. A reliable field-theoretic calculation of these matrix 
elements is the obstacle to a quantitative theory of hadronic weak decays. 

Using Fierz identities, the four-quark operators in the effective Hamiltonian may be 
rewritten in various forms. It is particularly convenient to rearrange them in such a way 
that the flavour quantum numbers of one of the quark currents match those of one of 
the hadrons in the final state of the considered decay process. In the case of B ^ D**7r 
transitions, omitting common factors, the various decay amplitudes may be written as 

^(50 ^ = (ci + ^){D**+7i-\{du){cb)\B'') 

+ 2c2{D**+n-\{dtau){ctab)\B'^), 
A{B'^^D**\^) = (c2 + —){D**\'^\{cu){db)\B'^) 

+ 2ci (Z}**V|(ct,n)(da)|5°), 
A{B^ D**\-) = A{B^ D**+7i-) - V2A{B^ D**V) , (11) 

where ta are the SU(3) colour matrices. The last relation follows from isospin symmetry 
of the strong interactions. The three classes of decays shown above are referred to as 
class-1, class-2, and class-3, respectively 0. 
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The class-1 amplitude A{B^ D**^7i ) contains the "factorizable contribution" 

Afact = {n-\idu)\0){D**+\icb)\B') , (12) 

which can be calculated in terms of the pion decay constant and the —>■ 
transition form factors. It also contains other, nonfactorizable contributions, which can 
be accounted for by introducing hadronic parameters ei and es such that 

{D**^n'\{du){cb)\B') = [l + ei{fi)] A^ct , 
{D**+n-\{dtau){ctMB'') = l^sifi) Afact . (13) 

Then the class-1 decay amplitude takes the form shown in (^), i.e. A{B^ D**^7t~) = 
aiAfact, with ||-ll| 



ai = (ci(/i) + ^) [1 + £i(/i)] + C2(/i) esifi) . (14) 

The hadronic parameter ai takes into account all contributions to the matrix elements 
and is thus fi independent. The scale dependence of the Wilson coefficients is exactly 
balanced by that of the hadronic parameters £j(/i). 

In the case of i? — D^*^7i: transitions, a similar discussion can be done for the class- 
2 amplitude, which contains the factorizable contribution {D^*^^\{cu)\0) {7i'^\{db)\B^) . 
Nonfactorizable contributions can be accounted for by introducing a parameter 02, which 
has a similar structure as ai except for an interchange of ci and C2. In the present 
case of B ^ D**tt transitions, however, things are more subtle, since the factorizable 
contribution to the class-2 decay amplitude vanishes. The reason is that the p-wave 
charm mesons Di and D2 do not couple to the vector or axial vector currents, i.e. 
{D**^\{cu)\ 0) = 0. Nevertheless, the matrix elements of the four-quark operators need 
not vanish. Taking into account that there is only a single helicity amplitude for the 
decays considered here, we choose to normalize these matrix elements to the factorized 
class-1 amplitude Afact introduced in ([12|) and define 

-V2{D**\''\{cu){db)\B') = 5i(/i) Afact, 
-V2{D*V\ictaU)idtab)\B') = i58(/i)^fact. (15) 

Then the class-2 amplitude takes the form —\^A{B^ D**^n^) = a2v4fact, where 

d2 = (c2(/i) + + ci(/^) Ssif^) ■ (16) 

Additional insight can be gained by combining these results with the l/N^ expansion 
0. At a scale /x = 0{mb), the large- A^'c counting rules of QCD implyQ ci = 1 + 0{1/N^) 



^For scales much lower than mt, the counting rules for the Wilson coefficients Ci(/i) are spoiled by 
large logarithms. 
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and C2 = 0{l/N() for the Wilson coefficients, and ei, 5i = 0{1/N^) and eg, = 0{1/Nc) 
for the hadronic parameters. Using these results, we find 



Hence, for the class- 1 amplitude we recover the result ai ~ 1, which also follows from the 
colour-transparency argument. The class-2 amplitude, on the other hand, is governed by 
a nontrivial hadronic parameter 6s{mf,) of order l/N^, which is process depend and will, 
in general, take different values for Di and ■ For the ratios of the various hadronic 
decay rates, we obtain 

r(5° ^ D**07r°) 1 2 



r(50 ^ D**+n-) 2 

T{B- D**°TT-) 
r(i?o ^ D**+n-) 



l + 6s{m,)\\ (18) 



Whereas the ffist ratio may be used to extract the magnitude (but not the phase) of 
^si'iTT'b), the second one determines the corrections to Bjorken's relation (||) for the case 
of class-3 decays: 



T{B- D 



dT{B- D**H-iy)/dr 



6n'f^\Vudr\l + 5s{m,)\\ (19) 



g2=0 



Note that the final state D**°7r is a pure isospin / = | state. Therefore, the class-3 
decays are unaffected by (elastic) final-state interactions. 



3 Semileptonic decay amplitudes for B — > D**i v 

The theoretical description of semileptonic decays involves the hadronic matrix elements 
of vector and axial vector currents between heavy meson states. These matrix elements 
can be parametrized as 



{D;{v',e)\c-f>'b\B{v)) = ikv{w)e 



(20) 



where w = v -v', and a mass-independent normalization of meson states is implied. The 
polarisation vector of the spin-one state Di satisfies e* ■ v' = 0; the symmetric, traceless 
Rarita-Schwinger spinor of the spin-two state D2 obeys the constraint e*^°f^ = 0. 
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The heavy mesons can be described by spin wave-functions with well-defined trans- 
formation properties under the Lorentz group and heavy-quark symmetry. To leading 
order in I/ttlq, expressions for the B — > D** transition matrix elements are readily 
obtained using the covariant trace formalism developed in Ref. ill]]. All form factors 



are proportional to a universal function r3(w), with calculable coefficients including the 
short- distance corrections to the decay amplitudes [|T2|. Explicitly, we find 



V2fvAw) 



V3 



w] 



1 

1 



ky{w) 



K-l)Cir|(w), 
3Ci + 2(w + 1)^2] r|H, 
(2 - w)Ci + 2{w + 1)C3] r3 (w) 

{w + 1)CIt.{w), 

(w; + l)Cf r|H, 

Ci T-3 (w) . 



(21) 



Here Ci and Cf are the Wilson coefficients appearing in the heavy-quark expansion of 
the weak currents (explicit expressions for these functions, which depend on w and the 
heavy-quark masses, can be found in Ref. |T2[)- 

It follows from these results that, in the heavy-quark limit, the matrix elements of the 
weak currents vanish at the zero- recoil point w = 1, refiecting the fact that the ground- 
state B meson is orthogonal to the p-wave D** states. The important observation made 
in Ref. |]13| was that the leading l/mg corrections at zero recoil can be calculated in a 
model-independent way in terms of the masses of charm-meson states. In this reference, 
rather complicated formulae have been derived where the various power corrections to 
the form factors are parametrized in terms of unknown, subleading universal functions. 
However, all model-independent information can be incorporated if the result for the 
form factor fvi{w) in (PT|) is modified according to 



V2fv,{w) = {w + l){w -1 + 26) Ci 
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1 



X 



2 



D 



M 



0.29, 



D 



(22) 



so that /vi(l) no longer vanishes. There are many other sources of l/mq corrections; 
however, they do not yield contributions at zero recoil. Since the physical range of w 
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values is restricted between 1 and ~ 1.3, it is likely that these other corrections will have 
a subdominant effect. We shall thus neglect them in our analysis. 
The differential semileptonic decay rates are given by 



dr{B^Di£-u) Gl\V,i,?M% 



dw 



dT{B Dl 



487r2 



V] 



dw 



727r3 



rl {w^ -1)2 {w- ra) fc^i + {w^ ~ 1) (^2/^^^ + A^^a) 



+ - (1 - 2wr2 + rl) 



1)4 



(23) 



where r\ = Md^/Mb ~ 0.459 and r2 = Md*/Mb ~ 0.466. The maximal values of w in 
the two cases are wq = {1 + r\)/2ri 1.32 and wq = {1 + r2)/2r2 ~ 1.31, respectively. 
When expressing the form factors through the universal function T-i{w) by means of the 
relations (0) and (H), it is convenient to introduce a new function t{w) ^ C,t^{w). 
The Wilson coefficient Ci is included in this definition since only the product of rs {w) 
with one of the short- distance coefficients is scheme-invariant. What remains in the 
expressions for the semileptonic rates are ratios of Wilson coefficients, which are scheme- 
independent and free of large logarithms. For all practical purposes, it is sufficient to 
evaluate these ratios at order as- 

Of particular interest for the further discussion are the differential decay rates at 
maximum recoil, corresponding to = 0. We find 



dT{B D** 



V) 



where 



dw 

K{D,) 
K{Dl) 



Gl\Vr,?M% (I 



wo 



7687r3 



(24) 



1 



+ 



1 — ri 
Cl 1 - r2 



1 + ri 
2ri 



C2 , C's 



0.95 + 1.85(5. 



r2^ + ^ 



0.90, 



Cl 2r2 

and the short- distance coefficients are evaluated at Wq. Numerically, we obtain 
dV{B Dii'p 



(25) 



dw 



dr{B D* 



dw 



0.160 (1 + 1.944(5)2 X [r(1.32)]2 ps"^ 
0.394 X [r(1.32)]2ps"S 
0.136 X [r(1.31)]2ps"^ 



(26) 
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To calculate the total semileptonic rates requires an ansatz for the form factor t{w). 
Since the accessible range of w values is small, it is sufficient to adopt a linear approx- 
imation with a slope parameter p^: t{w) = r(l) [1 — p^lw — 1) + . . .]. Higher-order 
terms can be partially taken into account by reexpressing the results for the decay rates 
obtained in linear approximation through values of the function t{w) at intermediate 
points. In that way, we find 

T{B^ Dii-u) ^ (0.0158 + 0.06945 + 0.116952) x [r(1.23)]2 ps"^ 

^ 0.0462 X [r(1.23)]2ps"S 
T{B D; ru) ^ 0.0207 X [r(1.21)]2 ps-1 . (27) 

The results for the differential and total semileptonic rates in ( PB] ) and (^) do not 
include l/mq corrections except those proportional to the quantity 6, which are kine- 
matically enhanced and specific for B ^ Di transitions. From the analysis of power 
corrections for B D^*H~u decays it is known that the remaining l/mq corrections 
tend to be spin independent and thus cancel, to a large extent, in ratios of decay rates 



1^ . Therefore, we expect that ratios of the decay rates estimated above are accurate to 
about 20%. 



4 Implications and conclusions 

We are now in a position to compare the theoretical predictions derived above with the 
available experimental data. Consider first the branching ratio for the decay B" —>■ 
D^l~v. Multiplying the theoretical prediction for the total semileptonic rate in (|27| ) by 
the lifetime tb- = 1.65 ps, we find 6(5" ^ D^ru) ^ 7.9% x [r(1.23)]2. This value 
is about 3 times larger than the result obtained in the strict heavy-quark limit, where 
the kinematically enhanced terms involving 6 are neglected. Under the assumption that 
B{D'( D*+Tx-) = 2/3, the experimental value in (|) implies 3(5" -> D^l'v) = 
(0.56 ± 0.15)%. Comparing this with the theoretical result, we find that r(1.23) = 
0.27 ± 0.04 is required in order to fit the data. This value is in good agreement with 
theoretical predictions obtained using relativistic quark models (r(1.23) ~ 0.3 |T^) or 



QCD sum rules (r(1.23) ~ 0.22 [15 ). We take this agreement as an indication that 
the terms involving the quantity 6 in (^B|) and (^^ do indeed capture the dominant 
corrections to the heavy-quark limit. 

Our next goal is to understand the ratio of the hadronic and semileptonic widths for 
B ^ Di transitions. The experimental data in (|1|) and (|^) imply 

^■- rg-^o;^^r ""^°°^- '''' 



independently of the Di D*tt branching ratios. Combining (0), (|26|) and (|27|) , we 
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obtain the theoretical prediction (using dg^ = 2MbMcii dw) 

dT{B- Dl(.-v)/<\w 



Ri 



ud 



MbMd, 



11 + Ssinib) 



WQ 



r{B- D^i-u) 



0.32 \l + 5s{mb)\ 



r(1.32)^ 
r(1.23)^ 



(29) 



which is rather insensitive to the value of 6. Obviously, there is no problem to account 
for the data. Making the conservative assumption = 1.5 ±0.5 for the slope parameter 
of the function t{w), we find Ri ~ (0.20 ± 0.06) |1 + 58(mfc)p, which is in good agree- 
ment with experiment provided the nonfactorizable corrections parametrized by 6s are 
moderate in size. Solving for these corrections we obtain |1 + 6s{mh) \ = 1.02 ± 0.25, in 
accordance with the fact that 6s = 0{1/Nc). 

Consider next the ratio of the two semileptonic rates for B decays into Di and D2 
mesons. The available experimental data for this ratio depend on some D** D*% 
branching ratios that have not yet been measured directly. We define 



h 



B(D? 



3.3 ±0.8, 



(30) 



where the quoted numerical value is obtained under the assumptions stated in the intro- 
duction. The CLEO data in (0) imply 



R2 



T{B- 



(0.16±0.19)/i = 0.53 ±0.63^ 



or i?2 < 1-48 (90% CL). From (|27|) , we obtain the theoretical prediction 

1.31 /rfl.2lV " 



R'? 



0.48. 



(31) 



(32) 



1 ± 4.395 ± 7.4052 \^r(1.23) / 

which is in good agreement with the data. We have again assumed = 1.5 ± 0.5 to 
estimate the form-factor difference, which is a small effect in the present case. Although 
the experimental errors in the value of R2 are large, to reproduce the data requires the 
presence of the terms involving the quantity 6 |p.3|| . 

Consider finally the ratio of the two hadronic rates for B decays into Di and D2 
mesons. The CLEO data in (p imply 



T(B- 



D 



*o 



vr 



T(B- 



(0.54±0.26)/i = 1.8 ± 1.0. 



(33) 



Combining (0) with the expressions for the semileptonic rates at maximum recoil given 
in (p6|), and assuming again ~ 1.5 to estimate the tiny form-factor difference between 
the two cases, we obtain the theoretical prediction 



R^ 



0.86 



1 ± 1.945)^ 



l + 6r'\m,) 



0.35 



l+6r^\mb) 



l±5r^^(m,) 



(34) 
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which is significantly lower than the data. It must be stressed, however, that the ex- 
perimental errors are large, and the discrepancy between theory and experiment is only 
about 1.5 standard deviations. Indeed, it would be very difficult to reconcile the central 
experimental value in (|33D with the theoretical expectation. The least understood aspect 
of the theoretical calculation is the question about the size of the nonfactorizable con- 
tributions to the hadronic decay widths. However, to reproduce the data would require 
that 



2.27 ±0.63. (35) 



Since the deviation of this ratio from unity is of order l/N^ in the large- A'^c limit, such a 
large value seems very unlikely. This is even more so as we have shown above that the 
nonfactorizable contributions to the B~ n~ decay rate are very small. Another 

possibility would be to blame the discrepancy between (^) and (0) on the D** — > D*7i 
branching ratios, which have not yet been measured directly. To reproduce the data 
would require that /i^0.65±0.31, which is much smaller than the commonly assumed 
value quoted in (pil|). Besides the fact that such a small value would imply exotic (or 
at least not well understood) physics in the strong decays of p-wave charm mesons, this 
possibility is essentially ruled out by the good agreement of the theoretical prediction for 
the ratio R2 with the experimental value for that ratio derived assuming the standard 
value for h. Finally, one may ask whether the theoretical prediction (0) could be 
spoiled by l/mq corrections not included in our estimate of the semileptonic rates in 
(^). However, as we have argued above, those corrections tend to cancel in ratios of 
decay rates and are expected not to exceed the level of 20%. 

To summarize, whereas there is in general good agreement between theoretical expec- 
tations and experimental data for semileptonic and hadronic B decays into final states 
containing a p-wave charm meson, the experimental central value of the branching ratio 
for the decay B~ — >• 1^2° cannot be accommodated by theory. We predict that future, 
more accurate measurements of this branching ratio will find a value 

B{B- Df -n-) ^ 0.35 B(fi- ^ L)° tt") ^ 4 x 10"^ , (36) 

which is about a factor 5 lower than the current central value reported by the CLEO 
Collaboration |jl]]. A similar conclusion has been reached recently in Ref. |l^ from a 
calculation based on the non-relativistic quark model. If, on the other hand, the current 
central value were conffimed, this would pose a serious problem for the theory of B 
decays. 

Acknowledgements: We are grateful to J. Bjorken, M. Beneke and B. Stech for helpful 
discussions. 
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